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ABSTRACT 



A theoretical model is derived to calculate the specific acoustic impedance of the 
absorptive material covering the walls of a cavity. This model will allow the exper- 
imental determination of the specific acoustic impedance from the measurement of the 
reverberation time in a water-lilled cavity. The model assumes a wall of low absorption. 
It can not be used for rigid or pressure release walls and grazzing incidence is excluded. 
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I. INTRODUCTION 



A. MOTIVATION 

The design of coating materials that reduce the reflection of water borne sound is 
an important problem in underwater acoustics. Since it is difficult to theoretically predict 
the complex reflection coefficient of such a material, the reflection properties must be 
determined experimentally. 

At high frequencies it is possible to measure the reflection coefficient of a slab of 
material by measuring the pressure after reflection of a narrow beam of sound incident 
at a given angle. However, this method becomes more difficult to use as the frequency 
of the sound is lowered. 

Determining the acoustic impedance of a material by measuring the standing wave 
pattern in an "impedance tube", while the standard procedure in air borne sound, is very 
difficult for water born sound, because of the impossibility of approaching "rigid condi- 
tions" for the tube walls. This is because the "pc" for water is approximately 5,000 times 
higher than that for air and is much closer to that of any wall material. 

Therefore it is attractive to investigate the possibility of determining the acoustic 
impedance of a matciial from measurements of the reverberation time in a water-filled 
cavity with walls lined with the material. 

B. OBJECTIVE 

To derive a theoretical model which will allow the experimental measurement of 
acoustical properties of absorption materials from measurements of reverberation time 
in a water filled cavity with wall lined with material. 

C. METHOD 

The model is based on the normal mode theory in cavities. 

It will be shown that the reverberation time in a cavity can be related to the 
temporal absorption coefficient of an acoustic field, f rom this fact, a relation between 
this coefficient and the complex specific acoustic impedance (resistance and reactance) 
will then obtained, under several assumptions and limitations. 

D. REMARKS 

No references were found using the proposed approach. In addition most of the 
published articles in this field are for air as the medium in the cavity. 
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II. LITERATURE REVIEW 



This is a summary of published articles which arc, in one way or the other, related 
to the subject. These summaries will give some idea about the accepted definitions and 
approximations relevant to this problem. 

A. CHAMBER FOR REVERBERANT ACOUSTIC POWER MEASUREMENTS IN 
AIR AND WATER [REF. 1] 

This was the only paper found in the Journal of the Acoustics Society of America 
where experiments made both in air and in water were reported, (in opinion corrob- 
orated by the authors). 

The motivation for this work was the interest in obtaining sound-power measure- 
ments directly, and the idea that it might be possible with the help of a reverberant 
acoustic chamber. It is stated that if the sound field consists of many randomly excited 
modes (diffuse field), the radiation impedance of the tank is not affected by the fluid re- 
action, because the randomly phased fluid pressures cancel. 

One of the major lessons to take out from this work is the fact that the reverberation 
time in the water, in a similar enclosure, is smaller than that of the air by a factor often. 
The author attributes this to the higher mean collision rate and higher sound speed in 
water. It is also noticeable that the air chamber gives a more uniform sound speed. This 
uniformity is important when the source radiation is concentrated in narrow frequency 
bands. 

The fluid loading factor is given by y = p 0 cj p P h2nf, where p„ and c 0 are the density 
and the sound speed in water, p„ and h the density and the thickness of the wall, and f 
the frequency. If y<l, the panels are very massive (hard walls) and specularly reflective; 
if y|>l, the walls are very' soft and specularly reflect sound as free surfaces. When y^l 
then the acoustic impedance of the wall is complex and the reflection is no longer 
specular, this the resistive and reactive parts of the wall impedance are function of the 
angle of incidence, this gives a more diffuse sound field. 

B. REVERBERATION TIME, ABSORPTION, AND IMPEDANCE [REF.2] 

Dowell presents a rigorous theoretical model to calculate reverberation time in a 
room in terms of the impedance of absorption materials on the wall. 



Ills approach is based on the coupling of the normal modes of a rigid wall room, 
through the damping due to wall absorption. 

The sound pressure in the room is given by the summation of the contribution of 
each normal mode p = p 0 c ; yP„F„IM n , where p 0 and c 0 arc the medium density and speed 
of sound, respectively, P n the coefficients in a normal mode expansion for pressure, F n 
for a rectangular room is co${n,nxlI. x ) co$(n y nyl L y ) cos^ipinz/ L z ) and M n = c x r. ) c l (for 
c„ — - 1 for n = 0 and c„ = 1/2 otherwise). l\ must satisfy the differential equation 

N 

P„+ co;P n + (A A p c .c 2 /l') ^P'CJAt, = 0 with A a being the area covered with absorption 

r—i 

materials, c o n the natural frequency of the mode and C„ = J ( F n FJz „) dA A jA A and the 
impedance of the material. 

From the assumption that the sound field is diffuse (i.e., uniform all over the room), 
the initial conditions for />. p, P n and P n can be calculated. In this case if we assume large 
impedance on the walls (zjp 0 c^> 1 ), all the coupling between modes are weak, and as 
P n = P„ = 0 at t = 0 for all modes except n = 0, the calculations can be made only con- 
sidering this mode. Other approximation can be made if we consider the initial potential 
and kinetic energy on each mode as being the same, then we can find the initial condi- 
tions for, respectively, P n and /'„. 

In the case of low absorption walls. i.e., large impedance, the damping coupling co- 
efficient is small (s n <;l), Dowell found the modal reverberation time 
T„ — 6 In 101 ’M n f ( FJ z„) dA A p 0 c] , except for 7'„, which value is 1/2 of T„. If a single mode 
dominates, the reverberation time is F n . 

Dowell claims that in the standard literature "the reverberation time is related to a 
random-absorption coefficient, the random absorption coefficient is related to the 
normal-absorption coefficient, and. finally, the normal-absorption coefficient is related 
to impedance" [Ref. 1: p. 1 S3], it does not consider room geometry, and considers uni- 
form absorption, lie considers his model as a generalization of the previous methods, 
since it computes reverberation time directly from wall impedance and geometric factors. 

C. SABINE REVERBERATION EQUATION AND SOUND POWER 
CALCULATIONS [REF.3] 

Various experiments show that the reverberation time changes with the amount and 
distribution of the absorption material in the cavity. This conflicts with the assumption 
commonly stated in the classical theory. 

Young says the measurement of the decay time instead of the reverberation time is 
easier and more direct. It is also advantageous when several absorption materials are 



tested on the same enclosure, because the results aic obtained by comparison and the 
use of the decay time eliminates the eventual effects of any bulk absorption in the me- 
dium. 

In order to obtain a diffuse sound field, required for the measurements, one must 
wait a while after the source is turned on, this is to give time for several reflections to 
occur. 

D. REVERBERATION TIME IN ENCLOSURES [REF.4] 

In this paper, the effect of the dependence of the sound absorption coefficient on the 
incident angle was tested. 

'I hey define a e as the dependence of the absorption coefficient a on the angle of in- 
cidence and it is given by a 9 = 1— | (Z/pc) cos O-lj(Zlpc) cos 0+1 1 : , where Z and pc are 
the acoustic impedance of the material and the medium, respectively, and 0 the angle 
of incidence. Furthermore the average absorption coefficient is given by 
a = f* i2 <x 0 cos Odd. 

This is a simulation of a two-dimensional enclosure, although its results permit some 
considerations about the physical problems connected to the levcrberation phenomenon. 
It supports the previous authors, who states that the reverberation time depends on the 
distribution of absorption materials. Even in perfectly diffuse fields, the reverberation 
time varies as much as 20% with respect to the Sabine's formula. 

Conclusion is that the angular-dependent sound absorption coefficient (a s increases 
as the incidence angle increases) has influence on the presence of specular reflections, 
what points to the believing that the sound absorption coefficient possesses the same 
property. 

E. OTHER LITERATURE 

All the other literature reviewed presented small or null interest for the problem. 
Most of them refer to work in air only and the approximation method do not consider 
normal mode theory and their application to experimental techniques is not convenient. 

From this literature the major points to take is the fact of the reverberation time of 
cavities with walls covered by absorption material be related to the acoustical impedance 
of the walls. All the other physical properties of interest, reflection coefficient and ab- 
sorption can be calculated from it. The article on the measurements of the reverberation 
time in air and water gives some feeling for the numerical results that one may expect 
from experiment. 
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All the mathematics of the theories was intentionally omitted from this shoit icview 
because it is of no interest for this work, all that we were looking for was the concepts 
and the mechanics of the problem. 



III. THEORETICAL MODEL 



A model for calculating the complex impedance of a low-absorption wall from 
measured reverberation times will be developed. It is assumed that the reader is familiar 
with normal modes in rectangular cavities; otherwise, section 9.7 of KFCS [Ref. 2: 
pp. 214-216] provides the necessary background to follow this development. All the sub- 
scripts (/, m, n) are dropped from the equations for simplicity, but their constant presence 
must not be forgotten. 



The geometry of the problem is presented in Fig. 1. 




A. REVERBERATION IN CAVITIES 

To measure the reverberation properties in a cavity, it needs to be completely Filled 
with sound, i.e., the sound field must be diffuse; the more normal modes excited, the 
better, so that a standing wave pattern is not identifiable. 

The number of normal modes N that are excited at a specific frequency f t is given 
by Knudsen [Ref. 3: p. 1 36]: 
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.'■-<■•(4)’ to 

Here V is the volume of the cavity and c the speed of sound in the medium. 

The equation governing the sound decay in this cavity with the sound turned off at 
t = 0 is 

P 2 =P( 0)e- ,l7 ‘ [2] 



where r, is the decay factor. 

The reverberation time (7) is defined as the time required for the sound pressure 
level to drop by 60 decibels (dB), so 



20 log = —60 



[3] 



or 



10 log — = 10 log e T,z ' = — 60 

^ no) 



which gives us, after some algebraic manipulation 

7 = 13.8t e 



[- 1 ] 



[5] 



B. DAMPED NORMAL MODES 

The lossy wave equation is 




-,2 

i c n 

c cT 



[6] 



where p is the complex acoustic pressure and c the complex speed of sound in the me- 
dium, which turn out to be real since we will ignore losses on the medium (relaxation 
time negligible). 

The Euler equation is 



= p 0 



dti 



ci 



where p 0 is the density of the medium and u the particle velocity. 



[ 7 ] 



n 

i 



As a solution to (6) try 



p = P cos (k^x+ipx) cos (ky.V+4> y ) cos [8] 

O) 

where t = ~ and the components of the complex wave number are: 



= k x +ja x 


[9] 


= k y +ja y 


[10] 


= k 2 -\-jx 2 


[11] 



where k x k y k, are the propagation terms and a x a y a, the absorption terms resulting from 
the wall expressed as a bulk absorption. T he complex frequency is written 



(0 = (O+jfi 



where p is the temporal absorption coefficient. 
Substituting (8) into (6) gives 




[ 12 ] 



[13] 



or equating real parts 






[14] 



and imaginary parts 



of} 

a x k x +or. y k Y +a 2 k 2 = — 

c 



[15] 



The particle velocity can be found by substituting the value ofp , from (8) into (7). 
Its components are: 



p k • [ 

u x sin (A r r+<p x ) cos (k Ji y+(p y ) cos(k ^+<p 2 )e , ~ 



[16] 



P fCy . 

M y =j — —cos(k x x+<l> x ) sh\{kyy+<p y ) cos (k^+cp.)/^- 



[17] 
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[18] 



u. =j-frTT co^(A x r+r/> x ) cos {kyV+fby) sin {kj+fije 1 - 1 

C. APPLICATION 

Assume that there is no bulk absorption and the wall at .v = l. x is the onh lossy and 
non-rigid wall. This implies </>,= 0 and a, y = a, = 0. Equation (15) simplifies to 

c a/? 

Mx=~r [i9] 

c 

Assume that the physical property which determines the absorption characteristics 
of the wall is the normal specific acoustic impedance, (there are no shear waves on the 
walls) and let this be defined as 



Iwall = (y X +j f l X )Po c 



[ 20 ] 



where y x and >] x are normal specific acoustic resistance and reactance coefficient. 

The normal specific acoustic impedance at the wall (say at jc = L x ) is by definition, 
the ratio of the acoustic pressure (p) to the normal component of the particle velocity 
vector («,), both evaluated at the wall. 



y r 



[ 21 ] 



At the wall, the wave and the wall impedances must match. This gives us the desired 
boundary condition 



r 



X- i , 



or, substituting £ (8), u x (16) and z wall { 20) 

jw~P 



Po I 



K x -p- x 



ctn {k x L x ) = {y x +jn x )Po c 



[ 22 ] 



[23] 



or rearranging. 



tan[(/< JC +y'a JC )L JC ] =j —j— 
aT j'Ix 



(0_ 

c 




k x +ja x 



[ 24 ] 
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we have the fundamental relation for the physical quantities important to this problem. 
The trigonometric identities 



sin(<7 + b) = sin a cos b+ cos a sin b 



cos(a+b) = cos a cos b— sin a sin b 



and defining 0 o — k,L, and <> = a ,L X , casts equation (22) into the form 

sin e 0 cos(/<5)+ cos 6 0 sin (J5) 

ta n{U 0 +jd) - cQs cos{jS) _ $in s - m( j 5) 

If <5<^l (i.e., low absorption) use of (19), and manipulation gives 



J 



(O 



k x c (Vx+j’h) 



1 JP «* 

'Mvr-x 



1 +jS ctn fl 0 
ctn 0 o -j <5 



[25] 



[26] 



[27] 



[28] 



Working the right hand side separately, multiplying it by ( tan 8J tan 0 o ) and by the 
complex conjugate of the denominator 



1+/ <5 ctn tan0„(l — a) . <5(l+taiT0„) 



+J 



1— y'<5taiT0 o l+<5* tan z 0 o v l+<) 2 tan‘0 o 



[29] 



where, since <5<^1 the term ( 1 — o’) 1. 

For the left hand side of(2Sj multiply by the complex conjugate of the denominator 



J 



O) 



'+4 V 



k K c(y x +j >] x ) 
i i 2 , O Vx( to 

k x c (Vx+>lx) L V 






+J 



. , ( JL 



)} 



Now introduce the incidence angle (measured from the normal) t //, 

k x 

COS t //, = — 



[30] 



[31] 



From equation (19) and knowing that cu/c = k 

k fi 



k x C 



[32] 
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Substituting a„ in 



co 



j*, 

k. 



P a x P 2 , 
~~r = ~ tan ^ 



CO 



So the right end side of (28) becomes (using the relation k t = •— cos <//,) 



cos 



i 

Phl+vl) { 



{ft 2 \ / [i 2 

n.- y\{ — tan V/,- j+/ — tan »//, 



P 



[ 33 ] 



[34] 



and equation (29) becomes 



O) 



cos 



>'Ai(v.x+^) - 



-y*( 77 tan2| Ai 



( 7 ) 



COS tib'x+fll) - 



-x( 77 tan^i) 



/*+’M 77 lan ‘^' 



tan 0 (O ( 1 — (5 2 ) _ <5(1+ tan^) 



[35] 



l+<5 tan‘(L 



+J 



l+<5 2 tan 2 P„ 



I. SOLUTION (general) 

Equating real and imaginary parts of (35), considering <5<^l . and (30), we get two 
equations for two unknowns y, and >/„: 



I 



('/.x ' x (i) tan Vj) 



cos i l'i{v x +’1 x ) 



tan B r 



1 +<5 tan Og 



[30] 



and 



cos ipiivl+vl) 



( P i \ <5(1+ tan 2 0 o ) 

■ Olx-Vx 77 tan h ) = — — — — 
v 7 l+<5 tan 0 o 



[37] 



Equations (36) and (37) can be more useful if manipulated. First divide one by 
the other and cross multiply 



>lx 



<5(1 + tan 0 o )~ — tan 0 o tan V? 

1 ? j 

tan 9 0 +S — tan~i/q(l + tan 0 o ) -> 



[ 3 S] 



The second equation can be obtained by the subtraction between (36) and (37): 



(]+(5 : lan 2 d 0 ) >/ A (l — tan 2 ^,)— y JC ( 
cos \pi(yl+>il) [ tan 0 o -«5( 1+ tan\)] 




[ 39 ] 



Now define a couple of parameters in order to simplify the expressions (38) and 
(39) for further work (note that 1/ tan 0 o = ctn 0 o ) 



As one may notice, "q" is mainly a function of the rigidness of the wall through 
the 9„ and "s" a function of the incidence angle i//, . 

Remember that, so far, the only assumption made is "low absorption" of the 
wall (<5<g 1 ). Let us see what that say about /?/ cm Start with equations (19) and (31) 



which, dividing both sides by co, multiplying the right end side by LJL,. and remember- 
ing the definition of d = o. t L x , becomes 



Since cos i/q < 1 , k and L^S^l , is clear that /?/«<? 1 

From this result we may induce that, if we avoid large values of i/q, than s<?U , 
i.e., the angle of incidence can not be close to grazing. 

Fig. 2 shows that it is also acceptable to consider q as being of order 3 if we 
exclude extreme values of ( ctn 0„-f tan 0 C ) , i.e., 6 0 can not be close to 0 (quasi-rigid wall 
boundary condition) or to n (quasi-pressure release wall boundary condition). 



P 2 , 
s = — tan i fo 

q = <5( ctn 0 o + tan 0 o ) 



[40] 




[41] 
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It can also be shown that, with some exceptions, s < q , so that 

— tan 2 i// ; < <5( ctn 0 o + tan 6 0 ) 

Using (42) and substituting k by oo/c gives 



£ 

10 



3 c 
to L x 



COS I p[ 



which when substituted in (43) and dividing both sides by 3 gives 



< ( ctn 6 0 + tan O 0 ) 

Assuming that kL,p 1 which excludes low frequencies (already excluded if we wish lots 
of normal modes excited in order to get a diffuse sound field) and since 
( ctn 0„+ tan O c ) ^ 2 (see Figure 2), the relation is true if we avoid grazing incidence (al- 
ready excluded). 

The assumptions and restrictions, so far, are: 
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- low absorption at the wall (n<l) 

- exclude close to grazing incidence (>//, £ nj2) 

- avoid rigid or pressure release wall, (0 C £ 0 and 0 o £ nj2) 

- s < q4,\. 

Note that the validity of the last relation can be extended when the dimensions 
of the cavity or the frequency are increased. 

2. RESULTS 

In this sub-section, expressions will be found for the normal resistance and 
reactance of the wall, using the normal incidence case as a check point. Insertion of the 
parameters s and q into (38) and (39) yields 



y.x _ izl 

} lx 1 ~\~sq 



[46] 



'bc(l--s) - y*(l+s) = (vl+’lx) cos I fo tan 0 o {\ - q) [47] 



Since (5qU and tan 0 o is not large, the term (l+<5 2 tan’OJ^l. 

Applying the approximations s < q<^\, which implies (q— s)<^l and s< 7<?1, it can 
be seen that the denominator of (46) goes to 1 and that (l-fs)^(l— s)~(l— <jr)~l The ratio 
yj>l 0 is much less than 1. It also tells us that both the resistance and the reactance have 
the same sign. 

Equations (46) and (47) become: 



>’x 




>i x -v 0 -(y 2 o +> i 2 o ) cos '/'/ tan e o 

Taking the value for y x from (48) and substituting on (49) yields: 
»/xC 1 -(‘7“-s)] = [((J-sf+llitil cos 1 1> i tan & 0 

Dividing both sides by q x and noticing that (q— i)<|l , we have an expression for t] x : 



[4S] 

[49] 



[50] 



>lx = 



1 

cos \pi tan 6 0 



[51] 
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Now substitute q, from (51) into (47) anti get y. 



< 7—5 

' x cos i//; tan 0 o 

or, plugging back the values for the parameters 5 and <7 (39): 

f! tan* </', 



/x cos I //; 



( ctn 2 d„+l)- 



10 



cos 1 / ’/, tan 



[ 53 ] 



The angle of incidence varies from 0 to n/2, so cos <//, is always positive, is a 
function of the mode number and the "rigidness" of the walls, so it can take any value, 
i.e., tan 0 o can be either positive or negative. Under these arguments, looking to (51) and 
(52), and remembering that (q—s) > 0 it is easy to verify that if tan Q, > 0 , then both y x 
and arc positive and they arc negative if tan 0 C < 0. This is consistent with the expec- 
tations stated earlier in this section. 

Applying the relation k, — w/c cos 1 //, to B 0 yields 



L x o) 

0 o = COS ij/ l 



[54J 



and a, can be related to /? by (31) 



1 } 1 

C COS \pi 






The temporal absorption coefficient /i can be obtained from the reverberation 
time T by direct comparison between equations ( 2 ) and (9) (on its maximum): 

pl = p]me~ P ‘ = [56] 



which leads to 




and with r, and T related by equation (5), we have 




[57] 



[5S] 



15 



So all the variables in (53) and (51) arc known or can be found from measure- 
ments, which permits the calculation of y, and ?/„ of the wall. 

3. CHECKS 



for d 



a. Main assumption 

First we can check the assumption cKl, to do so we may solve (51) and (53) 



\ U 0 



[59] 



The denominator is always greater than 0 (except for grazing incidence) , and it is larger 

V 

than 1 because we assumed it when checking for s < q, since <^1 our initial assump- 
tion is verified. 

b. Normal incidence 

We can start with (3S) and (39), considering again that (l+<5 2 tan 2 0J=l and 
if /, ■ = 0. Without any other assumptions that 5<1 and 0 o £(/?+ 1 )-/2 or n~, equation (38) 
becomes 



y x <5( 1 + tan 0 o ) 
Bx tan 6 0 



[60] 



Substitution of y x in (39) yields 

= ’/.xD+H-?) 2 ] tan e 0 

dividing both sides by >/„ and noting that (I— </)~l we get 






1 

tan 6 0 



[61] 



[62] 



which, substituting in (60) gives 






<7 

tan B 0 



= 5 ( ctn^+l) 



This results arc the same as (52) and (53) if we set i j/, = 0. 



[ 63 ] 
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c. Power reflection coefficient 

The power reflection coefllcicnt (R„) is given in KI CS JRef. 2: p. 139| by 
equation (6.45) 



K. 



(>v COS tffj-ly+til cos 2 !//; 

{y x c ° s 'h+iy+>il coi ' 2} i'i 



[64] 



We want it to be close to one (low absorption), so, dividing everything by ( cos 2 ^,) 



(y x - sec <l'j)) 2 +t?l—(y x + sec ^)) 2 +>/$ 



[65] 



This tells us that t/, can take any value, and 

yl-2 y x sec i/q+ secV/-y.t+2y x sec >//,■+ sec 2 i//, 



[ 66 ] 



01 



-Vjt sec 'pj—y x sec >/q [67] 

which becomes 

eos i//, [68] 

From this relation we find that y, despite of being much smaller than >]„ has 
to be small. Recalling equation (46) one may see that 

q—s 

y x = tttt n x < cos >Pi [ 6<) ] 

I +sq 

which, after some manipulation, and taking into account the approximations earlier 
made, tells us that rj x has to be also much smaller than ( cos i/»,). So we expect, for low 
attenuation at the walls, low resistance and low reactance. 
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IV. CONCLUSIONS 



Using normal mode theory in cavities it is possible to calculate the normal acoustic 
resistance and reactance as functions of the reverberation time, the frequency of the 
acoustic field, the speed of sound in the medium, the dimensions of the cavity, and the 
incidence angle. 

To measure the reverberation time, a dilfuse field is required, i.e., a fairly uniform 
pressure is desired throughout the cavity. 1 lie reverberation time can be related to the 
temporal absorption coefficient and it is expected to increase with the increasing di- 
mensions of the cavity, the decreasing absorbtion of the walls, and the decreasing pc of 
the medium. 

Under the following limitations: 

- low absorptive wall 

- avoid rigid or pressure release walls 

- exclude angles of incidence close to grazing 

the results obtained as a function of determinable parameters are: from (53), the normal 
specific acoustic resistance is 



Vx = 



CCOS't/q 



I. r (i) 

ctn'l — - — COS l/q 1 + 1 



tarft/q 



(O 



cos </q tan 



v> 

C COS </q 



[70] 



from (51), the normal specific acoustic reactance is 







[71] 



It is noticeable that, under the mentioned limitations, the resistive part has to be 
small and much smaller than the reactive part. 

The exclusion of the "close to grazing" angles of incidence is assumed as not being 
very important because, for a diffuse field, there are many excited normal modes, so the 
modes excluded by this limitation are relatively few. 
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Despite the limitations on tins model, it is a generalization of the theon picsented 
in the textbooks [Refs. 5, 6|. 
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